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Abstract 

We present a proposal for testing the prediction of non-equilibrium quantum field 
theory below the Schwinger limit. The proposed experiments should be able to de- 
tect a measurable number of gamma rays resulting from the annihilation of pairs in 
the focal spot of two opposing high intensity laser beams. We discuss the dependence 
of the expected number of gamma rays with the laser parameters and compare with 
the estimated background level of gamma hits for realistic laser conditions. 

Key words: quantum electrodynamics, non-pertubative limit, pair production, 
high intensity lasers. 
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1 Introduction 



Current generation of liigli-repetition and high-intensity lasers as well as 4*^^ 
generation light sources have opened the possibility to experimentally test 
quantum electrodynamics (QED) at the level where the details of multiple 
order expansions in the field propagators can be verified with measurable ob- 
servables. This subcritical field limit is of extreme scientific interest since it 
will allow progress of the physics research into completely new realms with 
the generation of novel and unexplored states of matter: electron-positron 
plasmas and excited vacuum states. Moreover, measured deviations from the 
predicted QED processes could indicate correction from quantum gravity or 
Lorentz-violations. These experiments can provide the required benchmark for 
cosmological vacuum particle production between the Planck and the GUT 
(Grand Unified Theory) era. In such cases the external field provided by the 
laser is replaced by the interaction of a massive scalar field with the back- 
ground space-time, but the governing equations of non-equilibrium quantum 
field theory (NeqQFT) still hold in the same form. 

In this paper we will discuss the theoretical framework implemented to calcu- 
late pair production at the subcritical limit and we will solve the governing 
NeqQFT equations in some idealized, yet representative, case describing the 
interaction of two high intensity laser beams in vacuum. In the second part 
of this paper, we will investigate the application of these approaches to a 
realistic experimental setup. In particular, we will discuss the possibility of 
observing gamma rays from pair annihilation in the laser focal spot and we 
will compare this number with the estimated background level. We will use 
as examples for the experimental configuration both the Gemini laser system, 
recently commissioned at the Rutherford Appleton Laboratory (UK), and the 
current generation of petawatt lasers such as the Vulcan laser at the Ruther- 
ford Appleton Laboratory. We will also compare optical pair production versus 
x-ray Free Electron Laser (FEL) sources and discusses the differences among 
those approaches. 



2 Subcritical pair production 



Due to its intrinsic simplicity, vacuum pair production is the proposed exper- 
iment of choice to test QED in the subcritical field limit. It is well known 
from basic QED that high energy photon-photon scattering can result in the 
production of electron-positron pairs. This process and the corresponding re- 
verse interaction of electron-positron annihilation play an important role in 
determining the overall opacity of the interstellar and intergalactic medium 
which in turn relates to the correct estimates of the intrinsic luminosity of 



stellar objects. At the same time, pair production provides a mean by which 
large fluxes of positrons could be generated in the vicinity of active galactic 
nuclei. One of the most exciting results of 7-ray astronomy has indeed been 
the detection of the 0.511 MeV emission line from the Galactic center [T]. 
Massive compact and quasi stellar objects are also source of intense beams of 
optical and infrared radiation and high-order low-energy multiphoton inter- 
actions which result in pair production are also possible. Moreover, processes 
involving a massive neutral scalar field in a dynamical background (either due 
to an external semi-classical field or a space-time metric) are applicable to 
cosmological problems such as vacuum particle production at the Planck time 
or reheating after GUT scale inflationary expansion [2|3P] . Another process 
closely related, and described within the same NeqQFT framework, is the 
thermal radiation arising from particle production near the event horizon of a 
black hole, commonly known as Hawking effects 0, as well as the Unruh effect 
[S], which is seen in uniformly accelerated detectors at relativistic speed. Re- 
cently, suggestions have been brought forward that the Unruh effect could be 
detected with the currently available Gemini laser at the Rutherford Appleton 
Laboratory [7]. 

On the theoretical side, subcritical vacuum breakdown is non-perturbative. 
Solutions exist only for idealized configurations and experimental verification 
is important for the correct understanding of the process and its relevance to 
the total interstellar opacity. At the same time this work could provide the 
first high-density electron-positron plasma to test and simulate a variety of 
astrophysical environments. With the advent of chirped pulse amplification 
(CPA) techniques [8] and progresses in x-ray free electron lasers (FELs) |9] 
it now becomes possible to generate very large numbers of coherent photons 
(i.e., high electric fields) in both the optical and the x-ray wavelengths. For 
any astrophysical object we define ?compactness? as the ratio between the 
total heating divided by its physical size [10]. In our context we are often 
dealing with high compactness objects, where pairs are primarily created by 
photon-photon collisions, and the energy loss is negligible as all the pairs 
remain confined within the laser focal spot. For the production of an electron- 
positron pair, the center of mass energy of the two photons must exceed 2inc^ , 
which precludes the creation of the pair by the collision of two single optical 
or x-ray photons. In strong electromagnetic fields, however, the interaction is 
not limited to initial states with two single photons, but allows multiphoton 
processes [TT] 

Ni{hoJi) + N2{huj2) ^ e-^ + e- , (1) 



where Ni and N2 are large integers. Experimental verification of the collision 
between ~4 coherent optical photons with one gamma ray photon with en- 
ergy ~30 GeV (created by Compton backscatter of another optical photon 



against an high energy electron beam) and the corresponding production of 
electron positron pairs has been demonstrated at the SLAG facility [T2II3]. 
Here, instead, we are interested in the more extreme case of a vacuum break- 
down driven by large number of low energy photons A^i ~ A^2 = -^ 3> 1 and 
oji r^ U2 = u! <t^ m (we use, as customary, natural units where Ti = c = 1). 
This is an example of NeqQFT where quantum mean field approaches have 
been proposed |2|14|15] but need experimental validation. The basic of these 
approaches is the so-called quantum Vlasov equation. In spinor QED assuming 
an external semi-classical electric field, it is possible to show that, for fermions, 
the particle number operator satisfies an equation of the type PIT3] 
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which is known as a quantum Vlasov equation, and 

^l = m' + kl + (A;|| - eA)\ (3) 

with 

e\ = m' + kl, (4) 

4 = (A:,| - eA)\ (5) 

and k± {k\\) is the momentum perpendicular (parallel) to the linearly polarized 
electric field S = —A (in the Coulomb gauge). The total electron-positron 
number per unit volume is then obtained by integrating over all the momenta, 

^W = 2/(&AW. (6) 

The quantum Vlasov equation has a non-Markovian character given by the 
factor [1 — 2/k(t)] arising from quantum statistics as it takes into account the 
full history of the distribution function. It simply says that the pair production 
rate will be affected by the particles already present in the system. However, 
in the case of weak (subcritical) fields, such an effect can be often neglected 
[T7] . The non-Markovian character is also inherent in the phase oscillations 
represented by the cosine term. This is related to quantum coherence, resulting 
from the fact that when the two pairs are created, they are initially fully 
correlated {i.e., entangled). The time-scale for these quantum coherence effects 



to wash out is in the order of Tqu ~ 27r/r2k ~ 2-n jm [2][H] . In order for the 
statistical description of pair production to be vahd, this time must be shorter 
than the time required to produce the pairs. This, for small /c, can be estimated 
from Eq. ([2]) to be 



Tcl 



^k(t)ex(t) 
fik(t)e||(t) 



S- (^) 



In the semi-classical case we also need to assume that the external field remains 
approximately constant during particle generation, that is r^i < Tpi \T^ , where 
Tpi is the characteristic time associated to collective plasma fluctuations: Tpi = 
27r/ujpi = 2tt [m/e^navY^'^ with Hav the average pair density. 

Equation ([6]) is just a formal solution, and a few words are necessary in order to 
correctly interpret its meaning. This point has been discussed in the literature 
[2115] . and it stems from the fact that the number of pairs does not commute 
with the Hamiltonian (indeed it is not a constant of motion). This follows 
directly from the uncertainty principle. If we have N^p pairs, the uncertainty 
relation reads as 

AEAt = A(A/'epm)At ~ 1, (8) 



and the uncertainty in the particle number is AN^p ~ 1/mAt. This implies 
that the particle number is indeed a well defined quantity at asymptotic times 
(At —7- oo) or for very massive (classical) particles. On the other hand, this 
relation applies only for a system where particle production during the time 
interval under consideration is negligible. In the more general case, we need 
to assume that particles will be produced within the considered time interval. 
We thus obtain 
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At, (9) 



which, letting dN^p/dt ~ Nep/Tcu is minimized for 



The particle number is a well defined quantity only if the change in the num- 
ber of particles is small within the time we are considering. If not, we need 
to resort to higher order approximations. This renormalization technique is 
referred to as adiabatic regularization [2|14|16] . We note that AA^ep cannot be 
made arbitrarily small, and it is minimized for AN^p = (2/m^/^) \dNep/dt\ ' ~ 



2[Ni.pe£Y/'^ /m. In summary, the quantum Vlasov equation represents a phys- 
ical observable (the number of electron positron pairs) only if the hierarchy of 
times Tmi^Tci and rg„ < Td < Tpi is satisfied. 

We should note that the NeqQFT framework is not the only one been imple- 
mented in the calculation of subcritical pair production, and at present a large 
amount of theoretical work has appeared |18fl9ll2Uf21f22f23f24f25] . While the 
various approaches seem to converge for large electric fields [I?], some dis- 
crepancies in the predicted pair number are seen in the subcritical regime. 
On the other hand, recent work seems to demonstrate that despite theoretical 
techniques being very different, they are effectively equivalent solution of the 
same problem, with the differences arising only from the details of the numer- 
ical methods j26jj. Still remains the fact that the precise details of the vacuum 
breakdown mechanisms in full spatial and temporal resolution are not yet fully 
understood despite the pioneering work of Schwinger [27]. Techniques based 
on the worldline path integral [28|29] as well as calculation of the tunneling 
probabilities of virtual pairs from the Dirac sea [30|3T] have been successful 
in determining the pair production in simplified non-uniform field configura- 
tions. However, experiments in supercritical fields {i.e., in the Coulomb field 
of an ion) have shown contradicting results and the question is still open on 
whether the Dirac equation is applicable in these scenarios [25] — see also 
discussions on the Klein paradox [32] — and the necessity to use a multi-body 
second quantization formalism (as in the NeqQFT approaches) becomes clear. 
In simpler terms, as particles are created, their associated electric field adds 
to the external field, which then feeds back to the production of the next pair. 



3 Solution of the quantum Vlasov equation for idealized fields 



Several attempts have been made to solve the quantum Vlasov equation in 
cosmological regimes [2||3|H] . More recently, attention has been drawn to the 
fact that the new generation of laser and PEL facilities has now reached electric 
field intensities where the particle production could have observable effects. In 
the simplest case of a time invariant, spatially homogeneous electric field, the 
solution is well known. This was originally derived by Schwinger [2|27] and 
found that the pair number is exponentially suppressed: 

AAschwinger OC CXp ( -^T— J = exp (-'^-^j , (H) 



where Sc = m? /e = 1.3 x 10^*^ V/m is the Schwinger field. Since the criti- 
cal field corresponds to the electric field such that its work on two electron 
charges separated by a Compton wavelength equals their rest mass, to reach a 



sizeable rate of pair production we need to have £^£c- On the other hand, the 
subcritical field regime is defined by the weak field condition £ ^ Sc-, implying 
a negligible number of pairs being generated. Equation fllip is only valid for 
static fields. In case of dynamically variable electric fields, the pair production 
problem can be understood as a tunneling with an oscillating barrier. This 
enhances the probability of generation of pairs since the average barrier seen 
by the virtual pair is lower. 

A clear advantage of the quantum Vlasov approach is that it can be used 
to model the full temporal dependence of the particle number for any time 
t>Tci. A solution of Equation ([2]) for a sinusoidal, spatially homogeneous, laser 
field has been recently proposed [24J. In this paper, we will consider instead 
a different temporal dependence of a spatially uniform field at the laser focus 
of two counter propagating laser beams: 



S(t) = Sosmh^{ut), 



(12) 



for which the Dirac equation is exactly solvable and analytical approximations 
are easily obtained. If we assume that the pair production is modest, i.e., /k ^ 
1, and So <S Sc, then ^33j 
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Under the condition e£Q/mv <^ 1, which implies a semi-classical motion of the 
charges in the electric field, the pair number can be further simplified to [33] 
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Using the field flT2|) . the asymptotic (residual) pair number density becomes 
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where we have assumed u '^ m, and e± ~ fik (which is valid for weak fields). 
Moreover, we have taken fik ~ "^ for k<m and fik ~ ^ for k>m. We see that 
in this case, pair production is exponentially suppressed for subcritical fields. 
The exponential term q-'^'^™-/'^ is indeed equivalent to what obtained with other 



techniques [30f3T] . This confirms the fact that, even for oscillating fields, a 
significant number of pairs can persist only for fields close to the Schwinger 
limit. The specific functional form for the residual density is dependent on the 
exact time variation of the electric field, and different results are obtained for 
sinusoidal fields S{t) = Sosm{ut) [2^|33] . In the latter case, the residual density 
after one oscillation period is rir ~ {cEqvY /m?, which is again negligible in 
the subcritical regime. 



While the residual density is exponentially suppressed at asymptotic times for 
the idealized field of Eq. f ll2p . the pair density at finite times is significantly 
larger. The average pair density during the field excitation can be approxi- 
mated as 
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Differently from the residual number, the average pair density is not exponen- 
tially suppressed. Moreover, calculation assuming a sinusoidal field showed the 
same functional dependence apart from a numerical prefactor of order unity 
[2^|33] . This may indicate that the average pair density is not too sensitive on 



the details of the field fiuctuations. 

Until now we have considered spatially homogeneous fields. Real fields, how- 
ever, are not spatially uniform and variations are expected to occur on some 
macroscopic scale A. These effects are more easily estimated within the semi- 
classical tunneling probability calculation [30|3T] . Since in a spatially inho- 
mogenous field pairs are initially produced at the maximum of the field, if 
they move away from this point and the field drops too sharply, they may 
not gain enough energy to cross the barrier and become real particles. Thus, 
opposite to the case of time varying fields, spatial gradients tend to suppress 
pair production. It can be shown that in the subcritical regime this effect 
introduces a correction to the pair production number of the order J30]|3T] 
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where m/e£oA<l. 



4 Observable effects from pair production 



As we have discussed in tlie previous section, in the subcritical regime, pair 
production at asymptotic times is always exponentially suppressed, meaning 
that no residual pairs remains after the laser. On the other hand, there is a 
significant number of pairs during the time the electric field is switched on. 
Assuming that the laser has wavelength A, then the estimated total number 
of electron-positron pairs in the laser spot volume V ^ s'^X (where s> A is the 
laser spot diameter) is given by 
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where the scale of spatial inhomogeneities is given by the spot size (A ~ s). 

Such number of electron-positron pairs has a clear observable effect, namely 
the generation of gamma rays due to pair annihilation. If during the laser pulse 
the particle number is a well defined physical quantity, then collision between 
those particles are indeed possible. Since the pair number scales with the laser 
wavelength {i.e., the interaction volume at the focal spot), it shows that opti- 
cal lasers have some advantage over x-ray FELs. On the other hand, for very 
intense FELs we could have the opposite scenario where the spot volume is 
too small to generate a sizeable number of pairs during the evolution of the 
laser pulse, but the field intensity is large enough that a finite number of pairs 
remains at asymptotic times. Those pairs may lead to accumulation effects 
{i.e., interacting with the external field) and induce spontaneous pair produc- 
tion over several laser cycles [31] , which is a consequence of the non-Markovian 
character of the quantum Vlasov equation. While the quantum Vlasov equa- 
tion is a collisionless equation, if collisions are a small perturbation, then the 
momentum distribution of the pairs is unaltered and the resultant number 
of gamma rays is obtained by multiplying their distribution functions by the 
annihilation cross section integrating over all the momenta of the pairs [24] . 
However, since we want to explore the scaling of the laser parameters with the 
observable number of gamma rays, we will follow here an equivalent analytical 
approximation. The ratio of electron-positron collisions producing gamma ray 
annihilation is 
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(19) 



where a^ is the Thomson cross section, which applies at low energies compared 
to the rest mass, as in our case. The model is applicable if collective plasma 
phenomena take places on a scale shorter than the laser cycle, i.e., Tpi < ^n/u, 
thus allowing sufficient time for the particle to interact [T4f55] . This also shows 



that as the laser field decreases, annihilation processes becomes more probable. 
Physically, this means that pairs produced with smaller momenta are more 
likely to result in a collision event. The number of gamma rays emitted is thus 
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where we have used the fact that u ~ 2tt/X (which is exact for a sinusoidal 
wave). While the treatment presented here is far from being complete, the val- 
ues obtained with this approach are in agreement with the predicted number 
of gamma rays calculated by full integration of the pair distribution function 
from the quantum Vlasov equation [21]. Both approaches are, however, not 
self consistent, and a complete analysis will require the addition of a colli- 
sional sink term in the quantum Vlasov equation [T5|37] . which must then be 
coupled to the gamma ray production rate. Only in this way the full effects of 
entanglement and quantum statistics could be properly accounted for. 

In order to get a realistic value for the expected number of gamma rays, we 
also need to account for the fact that counter propagating beam geometries 
are experimentally difficult to realize (see, however Ref. ^38] for a suggested 
counter-propagating beam geometry). If 9 is the angle between the two beams, 
then this introduces a geometrical correction (1 — cos0)/2. Moreover, if the 
laser beam has a pulse duration tl, then gamma ray annihilation events will 
occur tlI'/2tt times during the laser shot. 

Bringing back the factors of c and h, the total number of expected 77 events 
during a laser pulse is then 
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with the electric field is expressed in terms of the laser intensity, Jq, as £q 
(2/iocJo)^/2 (in SI units). 



5 Photometric 



In this section we compare expected number of 7 photons from pair annihila- 
tion with respect to background noises. As shown in Tabled! we expect ~0.6 
annihilation events per laser shot, corresponding to ~10000-events in a 10 hr 



10 



experiment using the Astra Gemini laser available at the Rutherford Appleton 
Laboratory. In a full experimental week, this corresponds to 5 x 10"^ annihi- 
lation events producing two gamma ray photons. Coincidence measurements 
will be performed with high sensitivity large area Nal gamma ray detectors 
covering a solid angle of ~ 27r, with an absolute conversion efficiency > 0.08 
|39] . We can estimate a total detection of ~ 2 x 10^ events. In situ measure- 
ments to assess the background level within the laser area have observed 2060 
positron events in 10 hours, equivalent to 0.05 counts/s. Since the Nal detec- 
tors can be gated with integration time ~1 /is, the background level of cosmic 
ray hits can be minimized to ~0. Any sporadic background event could be 
further eliminated with a coupled anti-coincidence detector. 



The major source of noise in these experiments arises from bremsstrahlung 
photons emitted by electrons stripped from the residual gas in the laser focal 
spot. Since relativistic electrons will be produced at laser intensities Jq^IO^^ 
W/cm^, this corresponds to a much larger volume than the laser focal spot. 
For the Astra Gemini laser, at pressures ~ 10~^ mbar, we expect up to 10^ 
electrons being ejected by the residual atoms (mostly hydrocarbons and oxy- 
gen). If these electrons are all emitted in a narrow cone, the probability that 
each one of them collides with a residual atom before reaching the chamber 
walls (~1 m path length) is less than 10~^. During such a collision a gamma 
ray photon is emitted, corresponding to < 0.04 events detected per laser shot 
(0.002 counts/s). If the gamma detectors are all placed within 1 m from the 
laser interaction point and outside the stripped electrons path, no additional 
gamma ray event will be recorded, as electrons hitting the chamber walls will 
emit photons in the forward direction away from the detector units (as well 
as excluded by coincidence detection). 



We notice from Table 1 that the error in the number of pairs is substantially 
larger than 1. This implies that the (rest) energy of those pairs is undetermined 
with an error (for the Gemini laser) of 1.3 MeV. Alternatively, this results can 
be interpreted in the sense that only a fraction of pairs has materialized on the 
mass shell, but the rest are still virtual. However, since the detection efficiency 
of scintillators remains the same over the ~1 MeV range (centered at 0.511 
MeV) |10] , we would expect that at worst a count rate is reduced by a factor of 
2, to 0.025 counts/s (1030 positron events in 10 hours), but still significantly 
above background. 
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Table 1 

Operation parameters for current laser systems and expected 77 yield. A beam 
crossing angle of ^ = 135° has been assumed. The Astra Gemini and the Vulcan 
PW systems are both located at the Rutherford Appleton Laboratory. 





Astra Gemini 


Vulcan PW 


Wavelength (nm) 


800 


1064 


Pulse length (fs) 


30 


500 


Laser energy (J) 


15 


500 


Spot diameter (/um) 


5 


5 


Intensity (W/cm^) 


2.5 X 10^1 


5 X 10^1 


£0 (V/m) 


1.4 X IQi^ 


1.9 X 10^4 


Uav (cm"3) 


8.0 X 10^0 


1.6 X 10^1 


Nep 


1.6 X 10^° 


4.2 X 10^° 


Tmi (fs) 


9.9 X 10-1° 


5.1 X 10-1° 


Tqu (fs) 


8.1 X 10^6 


8.1 X 10-^ 


Tel (fs) 


1.2 X 10^2 


8.7 X 10-3 


AiVep 


2.6 X 10^ 


5.0 X 10^ 


2-KTpl/iy 


0.13 


0.22 


jV'i"'' 


0.63 


0.21 


Repetition rate 


every 20 sec 


every 1 hr 


N^ff after 10 hr 


10879 


805 



6 Conclusions 



We have presented a proposal to test subcritical pair production with high 
intensity lasers. Using the theoretical framework of NeqQFT we have shown 
that the residual pair density after the laser shot is exponentially suppressed, 
and the number of pairs remaining is negligible. However, for realistic laser 
conditions, there is a significant number of pairs during the field evolution 
and the observable effect of such pairs is the production of co-incident gamma 
rays. We have estimated for the Astra Gemini laser facility at the Rutherford 
Appleton Laboratory more than lO"! annihilation events during an experimen- 
tal day. Photometric analysis has shown that this number of events will be 
detectable with current instrumentation. We are proposing an experimental 
platform that could test, for the first time, NeqQFT models which are relevant 
to astrophysical and cosmological processes, and, at the same time, resolve is- 
sues with the current approximation schemes of non-perturbative QED. 
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